function [SE_ITT,SE_SNT]=power_ind(n,C,tau,sigma,pi,f)

%Initialize distribution statistics
mu_k=zeros(1,length(pi));
s_k=zeros(1,length(pi));

%Calculate distribution statistics
for i=1:length(pi);
mu_k(i)=pi(i) *f(i);
s_k(i)=(1-pi(i))*f(i);
end;

psi=0;if pi(1)==0; psi=f(1);end;

%Standard Errors
SE_ITT = zeros(1,length(pi));SE_SNT = SE_ITT;
for i=2:length(pi);
SE_ITT(i)=(n*tau*(1/psi+1/f(i))+sigma*(1/mu_k(i)+1/psi)).^0.5/(n*C)^0.5;
SE_SNT(i)=(n*tau*(1/psi+1/f(i))+sigma*(1/s_k(i)+1/psi)).^0.5/(n*C)^0.5;
end;

end
